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Abstract 

A model of fluctuations in the market price including many deterministic 
dealers, who predict their buying and selling prices from the latest price 
change, is developed. We show that price changes of the model is approxi- 
mated by ARCH(l) process. We conclude that predictions of dealers affected 
by the past price changes cause the fat tails of probability density function. 
We believe that this study bridges stochastic processes in econometrics with 
multi-agent simulation approaches. 
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1 Introduction 

Market prices have been analyzed by many researchers. It is a famous result 
that the probability density function of difference of the market prices follows 
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the power law distribution (Mantegna et al. 1995). On the other hand the 
stochastic processes of ARCH type are one of the most exciting studies in 
econometrics. Statistical properties of ARCH processes have been clarified 
by many researches. Namely a probability density of a dynamical variable of 
the ARCH process follows the power law distribution and its characteristics 
is applied to explain fluctuations in the real market (Engle 1982, Bollerslev 
1986, Nelson 1990). Recently economists and physicists are interested in 
multi-agent simulations of problems motivated economically (Takayasu et al. 
1992, Bak et al. 1997, Sato et al. 1998, Johnson et al. 1998, Lux et al. 
1999). This approach is to study an artificial market in which programming 
agents sell and buy their stocks and actually market price changes can be 
calculated by an interaction of many dealers. 

On the viewpoint of statistical mechanics we can consider investigating 
statistical properties of market prices, the stochastic processes of the ARCH 
type and multi-agent simulations to be macroscopic, mesoscopic and micro- 
scopic, respectively, (see fig. [[]). Our study will bridge the multi-agent 
simulation with stochastic processes of the ARCH type. 

In this article we develop a simple market model with many deterministic 
dealers. The dealers estimate their buying and selling prices from the latest 
change of the market price. We show that from market price changes of the 
model a stochastic process of the ARCH type can be derived. Let us give a 
brief outline of this article. In the second chapter we describe the process of 
ARCH type. In the third chapter we mention a simple model of fluctuations 
in market price based on dealers. In the fourth chapter we analytically derive 
an ARCH(l) process for price changes from the proposed model. The fifth 
chapter is devoted to the concluding remarks. 



The ARCH abbreviates autoregressive conditional heteroskedasticity, which 
has been introduced by Engle in econometrics (Engle 1982). The ARCH(q) 
process is formalized by 



2 The ARCH (q) processes 




■ z s 




(1) 
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where e s is an interesting variable, Z s is a stochastic variable, a s is called 
volatility. a; (i = . . . q) is a positive parameter. Consider g = 1 as the most 
simple case of the ARCH(q) stochastic process. It seems sufficient to discuss 
the case of q = 1 in order to see statistical properties of the ARCH process. 
Especially when q = 1 eq. flU) is written by 



— C"s ■ Z s 

a 2 s = a + aiel_ v 



(2) 



Eliminating a s into eq. (§) one can get an alternative expression, 



'ocq + aie^! • Z s . (3) 
Moreover eq. (|3]) is approximated as 



v/tto-^s (|e s _i| < -y/ao/ai) 

V^oT- |e a _i| ■ Z s (|e s _i| > Ja Q /ai) 



(4) 



This equation is equivalent to a random multiplicative process (Takayasu et 
al. 1997). Using the random multiplicative process theory (Sato et al. 2000) 
it is easy to prove that a probability density function of e s has power law 
tails, 

p(e) ex \e\-P-\ (5) 

where f3 is given by 

a? /2 <i^> = I- (6) 



3 Dealer model 

We show a brief explanation of the market model consisted of many simple 
deterministic dealers. As shown in fig. ^| this model can be separated into 
two parts; One is a market mechanism, which describes how to determine a 
market price from orders. The other is an algorithm of agents, which governs 
how to order a selling or buying on the market and how to modify their bid 
prices. In the following subsections we explain these parts, initial conditions 
and parameters and also results of numerical simulations. 



3 



3.1 Market mechanism 



Suppose that a market of a competitive buying and selling in which N dealers 
give their orders for a common order board. Si(t) and B^t) represent a 
selling price and buying price for the ith dealer at time t, respectively. It is 
assumed that buying prices and selling prices individually compete in this 
market. Namely the maximum buying price and the minimum selling price 
are effective in the market. Thus the condition for a trade to occur is given 
by the inequality 

L{t) = max{£(t)} - min{S(*)} > 0, (7) 

where max{B(t)} represents the maximum buying price in all the buyers, 
and min{S'(t)} the minimum selling price in all the sellers. Suppose that 
the market price P(t) is determined as an arithmetic mean of the selling 
price and the buying one when the transaction occurs. Otherwise the latest 
market price is maintained. Namely, 

_ f §(max{£} + min{S}) (L(t) > A) 
P{t) -\P(t- 1) {L[t) < A) • (8) 

3.2 Dealer algorithm 

We suppose that a seller goes on decreasing his expectation of selling price 
until he can sell a stock and that a buyer goes on increasing his expectation 
of buying price until he can buy a stock. A rule to modify his expecting price 
is given by 

Bi(t+l) = Bi(t) + \l + cAP P rev\ai{t), (9) 

where AP prev denotes the latest change of the market price. From the above 
assumption for the seller and the buyer the ith dealer is a seller when a,i(t) < 
and he is a buyer when a,i(t) > 0. A term of |1 + CiAP prev \ means that 
modification of his expectation of price depends on the latest market price 
change. Here q is a prediction coefficient dependent on the dealer. 

Suppose that all the dealers have a small asset. Thus it means that a 
seller/buyer changes his position into a buyer/seller after a trade. Because 
of the assumptions that each dealer keeps his position till a trading, an 
evolution rule of a(t) is given by 

/, , x [ —aAt) (a seller and a buyer when a trade occurs) „ x 
aAt + 1) = { * y ) , . , J . (10) 

n ; 1 Oi(t) (otherwise) v ; 
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3.3 Initial conditions and parameters 

-Bj(O), eij(O) and q are given by random numbers of which a range is [— A/2, A/2], 
[—a, a] and [— c*,c*], respectively. We put P(0) = and AP prev = 0. The 
dealer's rule is deterministic except initial conditions. This model has four 
parameters; the dealer number N, a for cij(t), difference between selling and 
buying price A and c* for q. 

3.4 Numerical Simulation 

Fig. |3] show a typical example of time series of market price P(t) simulated 
numerically with price changes AP(i) = P(t) — P(t — 1). Because an event 
of a change of the market price occurs in Poisonian we define a jump Ap as 
a market price change on which a trade occurs. Fig. |] shows probability 
density functions (PDF) of the jumps and their corresponding cumulative 
distribution functions (CDF), which is defined by 



From fig. |] we can find linear part meaning that CDFs follow power law 
distributions of the following form 



where (3 is a power law exponent. Moreover the power law exponent obviously 
depends on a value of c*. 

4 From Microscopic to mesoscopic 

As shown in fig. [5] Ap s represents a change of the market price on the sth 
trade, and n s a time interval between the sth trade and the s + 1st. We 
define M s as the buying price at the sth trade and m s as the selling price, 
respectively. From the assumptions of eq. (^) Ap s can be estimated as 




(11) 



P(> \x\) oc \x 



(12) 



= -(m s + M s + A)--(m s _ 1 + M s _ 1 + A) 

1 1 

= -(m s -m s _i) + -(M s -M s _i). 



(13) 
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The first term and the second term on the second line are approximated by 

m s -m s ^i = |1 + CiAp 8 -x\a>in> 8 -i, (14) 
M s -M s _i = |1 + CjAp a _i|a 3 -n s _i, (15) 

where the subscript i denotes the dealer who gives the lowest selling price, 
and j the dealer who gives the highest buying one. Substituting eqs. (|14]) 
and (^) into eq. ( |T3"D yields 

Ap s = hi + CiAps-xlaiUs-x + hi + CjAps-tlajUs-i. (16) 

When the dealer number N is large it can be assumed that dealers' prices 
of expectation are distributed uniformly. By taking square of eq. (|16"D for all 
the dealers and averaging over ensemble under the condition that Ap s _i is 
realized, we get the following equation, 

M) = ^« 2 >(l + (c 2 )ApU)(nU)- (17) 

From the assumptions of c, and Oj we have (c) = 0, (c 2 ) = c * 2 /3, (a) = 
and (a 2 ) = a 2 / 3. Assuming that (Ap) = 0, eq. (|T7j) becomes, 

(A^) = « 2 (l + ^Ap 2 _i)(^i)- (18) 

Eq. ([i8|) shows that the variance of Ap s relates to a realized value on 
s — 1 and it has the same form with eq. (0) without the term of (n 2 _i). 
Supposing that a' s is a stochastic variable with a zero-average and a normal 
variance we can also rewrite eq. ([l8|) as 



c* 2 



A Ps = yi + — Ap 2 _i« \/(n 2 -xWs- (19) 

Eq. (|19D becomes identical to eq. (0) with the relations a = 1, «i = c* 2 /3 
and o" s = ay ( n i-i) cr 's- 
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5 Conclusions 



We introduced the model of an artificial market with many deterministic 
dealers and showed the outline of the stochastic process of ARCH type. The 
ARCH(l) process can be derived from the changes of the market price in 
this model. We conclude that the fluctuation of the market price is approx- 
imated by the ARCH type process if each dealer changes his expectation of 
a buying/selling price proportional to the lastest market price change. We 
expect that our approach will bridge the stochastic processes of the ARCH 
type in econometrics with dynamical market models consisted of dealers, and 
our understanding about the basic properties of markets will be deepened. 
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Figure 1: The conceptual illustration of correspondence of scope to study 
methods. 
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Figure 2: The conceptual illustration of dealer model. The inputs of the 
market are orders from dealers. The output of the market is a market price. 
The input of an agent is a sequence of latest price changes. The output of 
an agent is a selling or a buying. 
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Figure 3: Time series of artificial prices (right). Time series of changes 
of the artificial prices (left). X axis represents a step t. Parameters are 
N = 100,O! = 0.01,A = 1.0 and c* = 390.0. 
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Figure 4: Semi-log plots of the probability density functions of the changes 
(right). Log- log plots of the cumulative distributions of the changes (left). 
Parameters are N = 100, a = 0.01, A = 1.0. 
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P(Ks+l))«- 




C(s-l) x (s) x (s+l) 



Figure 5: A conceptual illustration of a temporal development of the market 
price P(t). t(s) represents a step for the sth trade to occur. Ap s is a price 
change on that step. n s denotes steps between the s — 1th trade and the sth. 
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